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The parton distributions functions (PDFs) derived from the NNLO QCD 
analysis of existing light-targets deep-inelastic-scattering data are pre- 
sented. The NLO and NNLO PDFs are compared in order to analyze 
perturbative stability of the analysis and estimate impact of the higher- 
order QCD corrections. The main theoretical uncertainties and experi- 
mental uncertainties in PDFs due to all sources of experimental errors in 
data are estimated and used to assess corresponding uncertainties in the 
cross sections of other hadronic processes. 

A study of the hard processes in the hadron collisions requires knowledge of 
the parton distribution functions (PDFs). The PDFs determine the normaliza- 
tion of the cross sections and therefore very often the uncertainty in the PDFs 
puts limit on the precision of measurements. For this reason the estimation of 
the PDFs uncertainties and improvement of our knowledge of the PDFs was 
the scope of intensive studies in recent years [1] . Estimation of the different 
uncertainties in PDFs shows that the uncertainties in the ansatz of analysis 
(theoretical ones) dominate over the uncertainties propagated from the experi- 
mental errors in data used to constrain the PDFs (experimental ones) for wide 
kinematics [2] . One of the potentially dangerous theoretical uncertainty in the 
PDFs extracted in the NLO QCD comes from the higher orders (NNLO) cor- 
rections. The NNLO corrections are unknown for the most of hard processes, 
but have been almost completely known for the DIS with calculation of the 
splitting function moments up to the 12-th [3]. The remaining uncertainty in 
the DIS structure functions due to missing terms in the NNLO correction was 
estimated at the level of few percent for the realistic kinematics [4] . 

The approximate NNLO corrections of Ref. [4] were applied in the QCD 
analysis of the existing charged-leptons DIS data including data from the SLAC, 
BCDMS, NMC, HI, and ZEUS experiments [5]. The NNLO and the NLO gluon 
distributions obtained in this analysis (A02) are compared in Fig.l. Variation 
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Figure 1: The lc experimental error bands for the A02 gluon distributions 
obtained in the NNLO (solid lines) and the NLO (dashes) at different values 
of Q. The NNLO gluon distribution for Set 1 of MRST2001 at Q 2 = 2 GeV 2 
is given for comparison (dots). 



of the gluon distribution under account of the NNLO correction is generally 
within the experimental errors that indicates perturbative stability of NNLO 
PDFs. In the same plot we give the NNLO gluon distribution extracted from 
the MRST global fit of Rcf.[6]. At low Q the MRST gluon distribution is 
well below the A02 one and gets negative at x <~ 0.0001. The origin of this 
difference has not been clarified entirely. One of the possible reasons is certain 
inconsistency of the NNLO MRST fit due to the NNLO corrections have been 
applied for the DIS only, while the MRST fit includes the Drell-Yan (DY) and 
the jet production data as well. This inconsistency was explicitly demonstrated 
in study of impact of the NNLO correction on the DY cross sections [7] : The 
NLO DY cross sections calculated using MRST PDFs are in good agreement 
with experimental data, while the NNLO predictions are well above the data. 
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Note for the jet data used in the MRST analysis the NNLO corrections are 
presumably even more important than for the DY process. Therefore, for the 
NNLO PDFs extracted from the DIS data the higher-order QCD theoretical 
errors are under much better control than for PDFs from the global fits. 
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Figure 2: The measurements of FL by the HI collaboration [8] compared to the 
NNLO QCD predictions (solid lines: A02(lcading twist+high twist); dashes: 
A02(lcading twist); dots: MRST(leading twist)). 



The NNLO QCD predictions for the structure function Fl based on the 
different sets of PDFs are compared to the HI measurements of Ref.[8] in 
Fig. 2. The MRST predictions are everywhere below the data. Since for the 
kinematics of Fig. 2 Fl is defined mainly by the gluon distribution this might 
mean that the MRST underscore the gluon distribution at small x and Q. For 
the A02 prediction agreement with the data is much better. At x <~ 0.0001 
the agreement is achieved mainly due to the leading-twist (LT) term, which 
was calculated using the NNLO PDFs. For x <~ 0.00001 the main contribu- 
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Figure 3: The NNLO VF/Z production rates in the pp collisions at y/s — 
1.96 TeV compared to the preliminary results for Run II [10] . The area between 
dashes gives la band uncertainty in the calculations due to theoretical and 
experimental uncertainties in the PDFs. 



tion comes from the high- twist (HT) term in Fl, which was fitted in the A02 
analysis simultaneously with the leading-twist term. This HT term might be 
manifestation of a new phenomena expected at small x or can be caused by 
the conventional effects like missing resummation of large logs or, eventually, 
by a disagreement of data from separate experiments used in the combined fit. 
In any case the unusual behavior of F L at small x and Q must stimulate more 
precise measurements in this region. 

With a good control of theoretical errors the NNLO PDFs extracted from 
the inclusive DIS data provide good tool for the precision measurements on the 
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hadron colliders. The NNLO predictions for the W/Z cross sections obtained 
using the code of Ref. [9] with the A02 PDFs are compared to the experimental 
data in Fig. 3. The prediction error bands given in Fig. 3 includes both theo- 
retical and experimental errors in PDFs. The accuracy of predictions is much 
better than the errors in data, which are mostly due to uncertainty in the lu- 
minosity. Therefore these prediction combined with the measurements of the 
W/Z rates can be used for better monitoring of the rates of other processes 
including exotic ones. 
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